Neutrino signals from the formation of black hole: 
A probe of equation of state of dense matter 
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The gravitational collapse of a non-rotating, black-hole-forming massive star is studied by v- 
radiation-hydrodynamical simulations for two different sets of realistic equation of state of dense 
matter. We show that the event will produce as many neutrinos as the ordinary supernova, but 
with distinctive characteristics in luminosities and spectra that will be an unmistakable indication 
of black hole formation. More importantly, the neutrino signals are quite sensitive to the difference 
of equation of state and can be used as a useful probe into the properties of dense matter. The event 
will be unique in that they will be shining only by neutrinos (and, possibly, gravitational waves) 
but not by photons, and hence they should be an important target of neutrino astronomy. 



Massive stars of >10 solar masses (M*©) end their lives 
when the iron core is formed at the center 0. For the 
mass range of ~10-20M Q , the spectacular optical dis- 
play known as the supernova explosion occurs following 
the gravitational collapse of iron core and the launch of 
shock wave by the bounce of core at high densities. The 
successful explosion produces a proto-neutron star that 
emits a bunch of neutrinos for ~20 s during the forma- 
tion and cool down |2(, which was vindicated in SN 1987A 

m 

Stars more massive than ^20M© may have different 
fates. They have larger iron cores and will be intrinsi- 
cally too massive to produce a neutron star via the su- 
pernova explosion. Then, the outcome will be a black 
hole. The gravitational collapse of these massive stars 
is currently attracting great interest. This is mainly be- 
cause they are supposed to be associated with gamma 
ray bursts with long durations, one of the most energetic 
explosions in the universe. Although the central engine 
and possible simultaneous production of hyper-energetic 
Type Ic supernova (sometimes referred to as hypernova 
in the literature) are remaining to be a mystery, many 
researchers believe that the collapse of rapidly rotating 
massive stars and the subsequent formation of black hole 
are responsible for the phenomenon 0. 

Not all the massive stars may not be rotating so 
rapidly, though. According to recent analysis of light 
curves of supernovae and nucleosynthetic yields by 
Nomoto and his collaborators there appears to be 
a subset of supernovae with a progenitor mass of ^20- 
25M Q , which produce a markedly smaller amount of 56 Ni 
and, as a consequence, are substantially underluminous. 
They suggest that these faint supernovae are slow rota- 
tors. If this is the case, it is highly likely that there exist 
essentially non-rotating massive stars that cannot pro- 
duce even a dim optical display and form a black hole. 



They are the target of this Letter. 

Such a non-rotating, black-hole-forming collapse will 
be no less faint in neutrinos than the rotational coun- 
terpart. Following the core bounce and shock launch 
as in the ordinary supernova, there will be a phase of 
hyper-accretion through the stalled shock onto the proto- 
neutron star, lasting for ~1 s before black hole formation. 
As shown below, as many neutrinos as in the ordinary 
supernova will be emitted during this period with dis- 
tinctive characteristics in luminosities and spectra. Then 
the emission will be terminated rather abruptly after a 
black hole is formed. These features will be an unmistak- 
able indication of black hole formation. The event will 
be unique in that they will be bright predominantly by 
neutrinos. Hence they should be an important target of 
neutrino astronomy. 

In fact, this channel of black hole formation (failed su- 
pernova), in particular, its neutrino signal has not been 
studied in detail with hydrodynamical simulations so far. 
Most of the previous papers were concerned with qua- 
sistatic evolutions of mass-accreting proto-neutron stars 
3 or the so-called delayed collapse of proto-neutron stars 
a EJ El El [3 113 that is supposed to occur after a 
successful supernova explosion, triggered by some phase 
transitions during its cooling phase. (Note we do not 
consider supermassive or intermediate-mass black holes 
[idj in this Letter.) In some other papers 0,0, long 
term hydrodynamical simulations of failed explosion sim- 
ilar to ours were performed but their foci were different 
from ours. Hence this is the first serious and quantitative 
investigation into the phenomenon. 

The most important findings in this Letter are the fact 
that the neutrino signals from the gravitational collapse 
of non-rotating massive stars are highly sensitive to the 
difference of equation of state (EOS) and can be used as a 
useful probe into the properties of hot and dense hadronic 
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matter. EOS of dense matter is one of the most impor- 
tant ingredients in high energy astrophysics but is ar- 
guably difficult to obtain from first-principle calculations 
and, as a result, the theoretical predictions are subject 
to uncertainty. Hence many attempts have been made 
over the years to obtain some constraints on EOS. We 
propose in this Letter to add a new item to this attempt- 
list, that is, the detection of neutrino signals from the 
hyper-accretion phase preceding the black hole formation 
in the non-rotating collapse of massive stars. 

A number of such events is admittedly highly uncer- 
tain but might be a substantial fraction (20~40%) of 
that of the ordinary supernovae, depending on the mass 
range for the event and the initial mass function of mas- 
sive stars. Then this event will be as important a target 
for neutrino astronomy as the ordinary supernova is cur- 
rently thought to be. It is, therefore, a very urgent task to 
provide theoretical predictions for luminosities and spec- 
tra of neutrinos. 

- Models. - We adopt the presupernova model of 4OM 
by Woosley and Weaver [Tjj • The model contains the iron 
core of 1.98M Q , which is larger than the ordinary size of 
~1.4M Q and warrants the black hole formation without 
explosion as its fate. We use the profile of central part 
of this model up to 3.0Mq in baryon mass coordinate to 
describe the accretion of envelope matter for a long time. 

We follow the dynamical evolutions by a general rela- 
tivistic z^-radiation-hydrodynamical code that solves the 
Boltzmann equation for neutrinos (v e , P e , tW T and v^/ T ) 
together with lagrangian hydrodynamics under spherical 
symmetry 0, 0, |2jJ. The general relativity is espe- 
cially crucial in the current study, since the re-collapse 
is triggered by general relativity. We assume the spher- 
ical symmetry, targeting the massive-end objects of the 
faint supernova branch, which are supposed to be rotat- 
ing slowly. 

In order to assess the influence of EOS, we employ two 
sets of realistic EOS by Lattimer and Swesty (LS-EOS) 
[13 and by Shen et al. (SH-EOS) [13 They have been 
used in most of recent studies of the ordinary supernovae 
|20| and are the current standard in the society. The 
subroutine of LS-EOS is based on the compressible liquid 
drop model for nuclei immersed in dripped nucleons. The 
data table of SH-EOS is relatively new and constructed in 
the relativistic nuclear many body frameworks [2^ | with 
the nuclear interactions constrained by the properties of 
unstable nuclei |2j| that are available recently in radioac- 
tive nuclear beams facilities. The two sets of EOS are, as 
a result, different from each other in stiffness at high den- 
sities and symmetry energy in neutron-rich environment 
|2*fij . For example, the incompressibilities are 180 MeV 
and 281 MeV and the maximum masses of cold neutron 
star are 1.8M and 2.2M Q for LS-EOS and SH-EOS, re- 
spectively. 

It should be noted here that both sets of EOS take 
into account only the nucleonic degrees of freedom and 
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FIG. 1: Position of shock wave as a function of time after 
bounce in models LS (thin) and SH (thick). 

ignore possible existence of other constituents such as 
hyperons, pion- and kaon-condensations and deconfined 
quarks. We do not intend to justify the negligence (see 
e.g. H H [H [H E El for their importance in the 
proto-neutron star cooling). It is, however, mentioned 
that this is essentially the first attempt to explore the 
EOS-dependence of the neutrino signals from the non- 
rotating collapse of massive stars leading to the black hole 
formation. Considering also the inherent uncertainties in 
the theories of these constituents other than nucleons, we 
think it is reasonable to start with the most basic EOS 
with nucleons alone, which will also serve as a reference 
for the future work with other degrees of freedom taken 
into account. 

The weak interaction rates for neutrinos are imple- 
mented by following the formulation of Bruenn |26| . 
In addition included are the plasmon process and the 
nucleon-nucleon bremsstrahlung process which is calcu- 
lated consistently with the effective nucleon mass in 
the EOS table |2j|. The other recent developments 
of neutrino-matter interactions [ljj and electron-capture 
rates [3] are not included in the current computations 
to facilitate the comparison with the previous study on 
15M0 star (2jj. It has been demonstrated that the 
electron-captures are important to determine the mass 
of the bounce core and the improvement may give a dif- 
ference of ~ O.IMq [2^]. The implementation of these 
reaction rates and more consistent treatments of reac- 
tions with EOS will be an important issue for the future 
work. We think, however, the basic scenario from the 
core bounce and shock stall through the hyper- accretion 
onto the proto-neutron star and then to the final collapse 
to black hole would not be changed substantially. 

- Results. - We start the description of numerical re- 
sults by presenting the radial position of shock wave as 
a function of time after the core bounce (tpb) in Fig. ^ 
In both models, the launched shock wave reaches the 
maximum radius and then turns into recession around 
tpb ~100 ms owing to the ram pressure of falling mat- 
ter. The size of bounced core is similar in two models 
and does not give much difference in the early dynam- 
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FIG. 2: Average energies (upper) and luminosities (lower) of 
v e (solid), v e (dashed) and f M / T (dot-dashed) as a function of 
time {tpb) in two models LS (thin) and SH (thick). 



ics of shock wave. Clear differences appear in the reces- 
sion of shock wave and the shrinkage of central core after 
t p b=100 ms when the hyper-accretion phase sets in. It 
should be noted that the accretion rate for the present 
model of 40M Q (~1 M Q /s at t p ;,=0.4 s) is considerably 
higher than that (~0.2 M /s) for the canonical 15M Q 
model for supernovae. This fact results in much faster 
contraction of central cores in the former. 

In model LS, the shock wave recedes quickly down to 
~20 km. The central core contracts rapidly as its mass 
increases toward the maximum value for hot and lepton- 
rich configurations in stable equilibrium. At t p b=0.56 
s, a dynamical collapse finally sets in and the central 
core shrinks on a dynamical time scale. By this time, 
the enclosed baryon mass and gravitational mass inside 
the shock wave reach 2.10M Q and 1.99M , respectively. 
Within the next ~8 ms, the central core becomes com- 
pact enough to form an apparent horizon at ^5 km, 
which marks the formation of black hole. 

In model SH, on the other hand, the shock wave recedes 
rather slowly over ~1 s. The dynamical collapse starts 
when the enclosed baryon mass reaches 2.66M Q (gravi- 
tational mass 2.38Mq) at t p b=1.34 s. This remarkable 
difference in the durations of the hyper-accretion phase 
preceding the black hole formation is worth particular 
emphasis. It originates mainly from the difference in the 
maximum mass of the hot and lepton-rich core in stable 
equilibrium and, to lesser extent, from the difference in 
the accretion rates. Hence, if observed, this will provide 
us with invaluable information on the stiffness of EOS. 

This novel difference is most clearly reflected in the du- 
ration of neutrino emissions as demonstrated in Fig. [21 
where the average energies and luminosities of neutrinos 



are shown as a function of time {tpb). The end points 
in the figure correspond to the formations of apparent 
horizon, i.e. the births of black hole. Note, however, 
that the major decline of neutrino emission will occur 
a fraction of millisecond later when the neutrino sphere 
is swallowed by the horizon and will be recognized at 
the boundary (^6000 km) another ~20 ms later, when 
neutrinos outside the neutrino sphere have traversed the 
distance at the light velocity. Unfortunately, we cannot 
follow this termination of neutrino emission owing to nu- 
merical problems. We will have to implement a scheme 
to avoid both coordinate- and real singularities to handle 
this problem (see, e.g., However, it is stressed 

that this problem does not matter in this Letter. The 
point here is that the longer-term neutrino emissions dur- 
ing the hyper- accretion phase is more revealing. 

The time profile of luminosities right after bounce is 
similar to the ones in ordinary supernovae having the 
neutronization burst of v e and the rise of v e , v^/ r and 
Vfi/r- Luminosities afterward are dominated by the con- 
tributions from the accreted matter, which is heated up 
by the shock wave and further by compression onto the 
proto-neutron star surface. Since the accreted matter 
contains a lot of electrons and positrons, they annihi- 
late with each other to create pairs of neutrino and 
anti-neutrino of all species. They are also captured by 
nucleons to produce electron-type neutrinos and anti- 
neutrinos. This latter processes are responsible for the 
dominance of v e and v e as well as their similarity in the 
luminosity. 

The difference in the reactions also leads to the differ- 
ence in the radial positions of neutrino sphere and, hence, 
to the hierarchy of average energies shown in Fig. [21 The 
average energy of v^/ T and v^/ T is particularly a good in- 
dicator of the difference of temperatures in two models, 
having, for example, a higher average energy at t p b ~0.5 s 
in model LS owing to the faster contraction. It is remark- 
able that the luminosities and average energies increase 
by a factor of two or more toward the formation of black 
hole, which will be utilized for diagnosis of the present 
channel of black hole formation. 

The earlier contraction of proto-neutron star and for- 
mation of black hole in LS apparently comes from the 
softness of the EOS, which gives a smaller maximum 
mass of neutron star. In Fig. [31 we display the evolu- 
tion of central density as a function of time (t p b). Right 
after a small peak at bounce, the central density in model 
LS rises very quickly toward the final collapse in contrast 
to the much slower increase in model SH. Since the con- 
traction of proto-neutron star proceeds almost adiabat- 
ically, the temperature inside becomes enormously high 
owing to the compression. The peak temperature is at- 
tained off center near the inner core surface thanks to the 
shock heating. We show in Fig. [3 the time evolution of 
temperature at 0.6M Q in baryon mass coordinate, which 
is near the temperature peak. The temperature reaches 
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FIG. 3: Baryon mass density at center (left) and temperature 
at off-center (right) as a function of time (t p f,) in models LS 
(thin) and SH (thick). 

around 100 MeV at the beginning of re-collapse and ex- 
ceeds much over 100 MeV for both models while it rises 
much slower before re-collapse in model SH. 

It is to be noted that we have used LS-EOS and SH- 
EOS in very high density and temperature regimes that 
they may not have been originally meant for. (We have 
extended the SH-EOS table using the original framework 
for the current simulations.) This will certainly be an 
over-simplification. Exotic phases such as meson con- 
densations, hyperons or deconfined quarks are expected 
to appear at certain densities and temperatures. The 
current models will serve as a reference for comparison. 
If, for example, a new phase emerges at a certain point in 
evolution, the softening of EOS will cause immediately 
a dynamical collapse to black hole just like the ordinary 
delayed collapse. Hence the present simulations provide 
upper limits to the duration of neutrino emission before 
black hole formation. It should be emphasized again that 
the difference in EOS is already remarkable in the neu- 
trino signals without these exotic constituents, which is 
in good contrast to the proto-neutron star cooling. 

The current generation of neutrino detectors will de- 
tect thousands of neutrinos from the above-discussed 
events in our galaxy. Simultaneous observations of differ- 
ent neutrino flavors, from the initial burst of v e through 
the rise of luminosity and hardening of spectra up to the 
termination, at Super-Kamiokande, SNO and other facil- 
ities together with an appropriate consideration of neu- 
trino oscillation are indispensable to claim the black hole 
formation through the channel of current interest. The 
neutrino signal may be a complex folding of a couple of 
factors in addition to EOS. It is of urgent importance 
to provide more detailed theoretical predictions, partic- 
ularly the dependence not only on EOS but also on the 
mass and structure of progenitor, which we are currently 
doing systematically. 

In summary, we computed neutrino signals from the 
collapse of a non-rotating massive star of 40M© and 
found that they are remarkably sensitive to the differ- 



ence of EOS and, hence, can be utilized as a novel probe 
into the properties of hot and dense matter. The ter- 
restrial detection of such events in future will reveal the 
quantitative details of the accretion and contraction of 
proto-neutron star deep inside the star leading to black 
hole formation, and will give a new constraint on EOS at 
high density and temperature. 

The numerical simulations were performed at 
NAO/ADAC (wks06a, wknlOb) and JAERI, and were 
done partially at KEK (No. 126) and RIKEN. This 
work is partially supported by the Grants-in-Aid for 
the Scientific Research (14039210, 14079202, 15540243, 
15740160, 17540267) of the MOESSC of Japan, for 
the 21st century COE program "Holistic Research 
and Education Center for Physics of Self-organizing 
Systems", and for Academic Frontier Project of MEXT. 
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